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We studied nonequilibrium oscillations between left-handed electron neutrinos and 
nonthermalized sterile neutrinos. The exact kinetic equations for neutrinos, written 
in terms of neutrino density matrix in momentum space were analyzed. The evolu- 
tion of neutrino density matrix was numerically calculated. This allowed to study 
precisely the evolution of the neutrino number densities, energy spectrum distor- 
tion and the asymmetry between neutrinos and antineutrinos for each momentum 
mode. Both effects of distortion and asymmetry, which cannot be accounted for 
correctly when working in terms of particle densities and mean energies, and the 
depletion of electron neutrino state have been proved considerable for a certain 
range of oscillation parameters. The influence of nonequilibrium oscillations on 
primordial nucleosynthesis was calculated. Cosmologically excluded regions for 
oscillation parameters were obtained. 

1 Nonequilibrium neutrino oscillations 

We discuss nonequilibrium oscillations between weak interacting electron neu- 
trinos z/e and sterile neutrinos for the case when do not thermalize till 

2 MeV and oscillations become effective after i/^ decoupling. Oscillations of 
that type, but for the case of i/g themializing before or around 2 MeV have 
been already discussed in literatureErB We have provided a proper kinetic 
analysis of the neutrino evolution in terms of kinetic equations for the neu- 
trino density matrix in momentum space. The assumptions of the model are 
the following: (a) Singlet neutrinos decouple much earlier than the active neu- 
trinos do: > Tjf , therefore, in later epochs T^, < T^,^, due to the addi- 
tional heating of Vg in equilibrium in comparison with the already decoupled 
Vg. Hence, the number densities of Vs are considerably less than those of i/g, 
Ni,^ << N^^. (b)We consider oscillations between i/s and Ve, according to 
the Majorana&Dirac mixing schemeEl with mixing present just in the electron 
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sector Vi ~ Uuvi, I — e,s^ 



V2 — —sm{d)ve + cos(d)vs, 



(1) 



where Vs denotes the sterile electron antineutrino, vi and 1/2 are Major ana 
particles with masses correspondingly mi and TO2- (c)We assume that neutrino 
oscillations become effective after the decoupling of the active neutrinos, Tqsc ^ 
H for T < 2MeV. This puts constraint on the neutrino mass difference: Sm^ < 
1.3 X 10~^ eV^. (d) We require that sterile neutrinos have not thermalized till 
2 MeV when oscillations become effective. puts the following limit on 

the allowed range of oscillation parameters: lliltra sin{2'd)5m? < lO^^eV^. 

As far as for this model the rates of expansion of the Universe, neutrino 
oscillations and neutrino interactions with the medium may be comparable, 
we have used kinetic equations for neutrinos accounting simultaneously for the 
f npntrinns into expansion, oscillations and interactions with 



participation. 

the mediumliiltjflU We have analyzed the evolution of nonequilibrium oscillat- 
ing neutrinos by numerically integrating the kinetic equations for the density 
matrix in momentum space for the period after the decoupling of the elec- 
tron neutrino till the freezing of neutron-proton ratio (n/p-ratio), i.e. for 
2 MeV > T > 0.3 MeV. We considered both resonant 6m? = m^ — m\ < Q 
and nonresonant Sm^ > oscillations. 



2 Kinetics of nonequilibrium neutrino oscillations 

The kinetic equations for the density matrix of the nonequilibrium oscillating 
neutrinos in the primeval plasma of the Universe in the epoch previous to 
nucleosynthesis have the form: 

^ = Hp ^ + 1 [Ho, p{t)] + 1 [n„u,p(t)] + o (hL) , (2) 

where p is the momentum of neutrino and p is the density matrix of the massive 
Majorana neutrinos in momentum space. 

The first term in the right side of Eq.|^ describes the effect of expansion, the 
second is responsible for oscillations, the third accounts for forward neutrino 
scattering off the mediums Tio is the free neutrino Hamiltonian: 



Tltei transitions between different neutrino flavours were proved to have negligible 
effectsEl 
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while TLint ~ a V is the interaction Hamiltonian, where = U*JJje, V ~ 
Gp (±i — Q/My^,), and in the interaction basis has the form 

LR_f V 



^--^0 )■ 

The first 'local' term in V accounts for charged- and neutral-current tree-level 
interactions with medium protons, neutrons, electrons and positrons, neutrinos 
and antineutrinos. It is proportional to the fermion asymmetry of the plasma 
L = Lf, which is assumed of the order of the baryon one 

Lf i i- ^ = f3T^. 5 

The second 'nonlocal' term arises as an W/Z propagator effect^ Q ^ E^, T^. 
The two terms have different temperature dependence and an interesting in- 
terplay between them during the cooling of the Universe is observed. The last 
term in the Eq.H describes the weak interactions of neutrinos with the medium.| 
We have analyzed the evolution of the neutrino density matrix assumed 
that oscillations become noticeable after electron neutrinos decoupling. So, 
the neutrino kinetics down to 2 MeV does not differ from the standard case, 
i.e. electron neutrinos maintain their equilibrium distribution, while ster- 
ile neutrinos are absent. Then the last term in the kinetic equation can 
be neglected. The equation results into a set of coupled nonlinear integro- 
differential equations for the components of the density matrix:| We have 
numerically calculated the evolution of the neutrino density matrix for the 
temperature interval [0.3, 2.0] MeV. The oscillation parameters range studied 
is Sm^ G ±[10^^*^,10^^] eV^ and i9 £ [0, 7r/4]. The baryon asymmetry /3 was 
taken to be 3 x IQ-^". 



3 Nucleosynthesis with nonequiUbrium oscillating neutrinos 

We analyzed the influence of nonequilibrium oscillations on the primordial pro- 
duction oi^He. The effect, ofxisci|llations on nucleosynthesis has been discussed 
in numerous publicationsErl§ll3ll3 Here we provided a detail kinetic calculations 
of helium abundance for the case of nonequilibrium oscillations in medium. The 

For example, for the weak reactions of neutrinos with electrons and positrons e+e~ ^ 
ViUj, e^Vj — > it was explicitly written in Refs. 10,12. 

''For the case of vacuum neutrino oscillations these equations were analytically solved in 
Ref. 12. 
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kinetic equation describing the evolution of the neutron number density in mo- 
mentum space n„ for the case of oscillating neutrinos <-> Vs was numerically 
integrated for the temperature range of interest T e [0.3, 2.0] MeV . 

4 Results and conclusions 

Our numerical analysis showed that the nonequilibrium oscillations can con- 
siderably deplete the number densities of electron neutrinos (antineutrinos), 
distort their energy spectrum and produce neutrino- antineutrino asymmetry 
that may grow at the resonant transition and may change considerably the 
evolution of neutrino ensembles. The effects of nonequilibrium oscillations on 
nucleosynthesis may be considerable for certain range of oscillation parameters. 
The results of our study are as follows: 

(a) As far as oscillations become effective when the number densities of 
Ve are much greater than those of Vg^ N^^ 3> iV^^, the oscillations tend to 
reestablish the statistical equilibrium between different oscillating species, As 
a result N^^ decreases in comparison to its standard equilibrium valuel3"Ll 

The depletion of the electron neutrino number densities due to oscillations 
to sterile ones leads to an effective decrease in the weak processes rates, and 
thus to an increase of the freezing temperature of the n/p-ratio and corre- 
sponding overproduction of the primordially produced ^_ff e. 

(b) For the case of strongly nonequilibrium oscillations the distortion of 
the energy distribution of neutrinos may be considerable.'! The evolution of 
the distortion is the following: First the low energy part of the spectrum 
is distorted, and later on this distortion concerns neutrinos with higher and 
higher energies (Fig. [^) . 

This behavior is natural, as far as neutrino oscillations affect first low 
energy neutrinos, Fqsc ~ 6rn? / Ei,. The naive account of this effect by shifting 
the effective temperature and assuming the neutrino spectrum of equilibrium 
form gives wrong results for the case 5m? < 10~'^eV^. 

The effect of the distortion on primordially produced helium is as follows. 
An average decrease of the energy of active neutrinos leads to a decrease of the 
weak reactions rate, F^, ~ and subsequently to an increase in the freezing 
temperature and the produced helium. On the other hand, there exists an 
energy threshold for the reaction + p —>■ n + . And in case when, due 

"^This effect was discussed for a first time by Dolgov0 but unfortunately, as far as the 
case of flavour neutrino osciUations were considered and the energy distortion for that case 
was showf to be negUgible, it was not paid the necessary attention it deserved. In paper by 
Kirilovafj it was first shown that for the case of Ue ^ Us vacuum osciUations this effect is 
considerable. 
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Figure 1: The figure shows the energy distortion of active neutrinos x'^pi^i^^x), where x = 
El, IT , for the case of nonequilibrium neutrino oscillations, hm?" = — 10~**, i3 = n/S at 
different temperatures: T = 1 MeV (a), T = 0.7 MeV (b), and T = 0.5 MeV (c). 



to oscillations, the energy of the relatively greater part of neutrinos becomes 
smaller than that threshold the n/p- freezing ratio decreasesllll The numerical 
analysis showed that the latter effect is less noticeable compared with the 
previously described ones. 

(c) Other interesting effect revealed by our approach is the generation of 
asymmetry between Ve and their antiparticles. The possibility of an PfYm- 
metry generation due to oscillations was discussed in many papersEEtJEr EtJEl 
Our approach allowed precise description of the asymmetry and its evolution, 
as far as working with the self- consistent kinetic equations for neutrinos in 
momentum space enabled us to calculate the behavior of asymmetry at each 
momentum. The calculated result may differ considerably from the rough esti- 
mations made by working with neutrino mean energy and with the integrated 
quantities like particle densities and the energy densities. The asymmetry effect 
is noticeable only for the resonant case. Even when the asymmetry is assumed 
initially negligibly small (of the order of the baryon one), i.e. ~ 10^^'^, it may 
be considerably amplified at resonant transition due to different interactions of 
neutrinos and antineutrinos with the CP-odd medium. The value of the asym- 
metry may grow by several orders of magnitude, oscillating and sign changing. 
Even in case when the value of the asymmetry does not become considerable 
enough to have some direct noticeable effect, on primordial nucleosynthesis 
for example, the asymmetry term at the resonant transition determines the 
evolution of the neutrino density matrix. It effectively suppresses the resonant 
transitions of active neutrinos (antineutrinos) thus weakening neutrino deple- 
tion at resonance. For some model parameters this effect consists 20% of the 
previously discussed. 

The asymmetry calculations showed a slight predominance of neutrinos 
over antineutrinos, leading to decrease of helium. The greater effect of the 



5 



Xn(t) i\ 
0,3 




matter oscilation - 
vacuum oscilation - 
standard model 



2 4 6 8 t(scc) 

Figure 2: The curves represent the evolution of the neutron number density relative to 
nucleons X„{t) = Nn{t)/{Np + Nn) for the nucleosynthesis model with vacuum nonequilib- 
rium oscillations and for the case of nonequilibrium oscillations in medium, 5m^ = — 10~**, 
1? = 7r/8. For comparison the curve corresponding to the standard nucleosynthesis model is 

shown. 

asymmetry is, however, the change in the evolution of the neutrino ensem- 
bles during and after the resonance resulting to a relative increase of both 
neutrino and antineutrino particle densities. This may lead to a noticeable 
underproduction of helium (up to 10% relative decrease). 

The total effect of nonequilibrium neutrino oscillations is overproduction 
of helium in comparison to the standard value. The results of the numerical 
integration are illustrated on Fig. |^, in comparison with the vacuum case and 
the standard nucleosynthesis without oscillations. 

From numerical integration for different oscillation parameters we have 
obtained constant '^He contours. We have used the 4% relative increase in the 
primordially produced helium to obtain the exclusion region for the oscillation 
parameters (Fig. 0). 
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Figure 3: Exclusion regions for oscillation parameters are shown for the case of resonant 
5m? < and nonresonant 5m^ > neutrino oscillations. The curves correspond to helium 

: 0.245. 



For the cases when the energy distortion and asymmetry are considerable 
we have Oibtained an order of magnitude stronger constraints than the cited in 
literatureo Therefore, in conclusion we would like to stress once again, that 
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in case of nonequilibrium neutrino oscillations working with the exact kinetic 
equations for the density matrix of neutrinos in momentum space is necessary. 
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